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Abstract 

The surface structure of Bi(llO) has been investigated by low-energy electron diffraction (LEED) 
intensity analysis and by first-principles calculations. Diffraction patterns at a sample tempera- 
ture of 110 K and normal incidence reveal a bulk truncated (lxl) surface without indication of 
any structural reconstruction despite the presence of dangling bonds on the surface layer. Good 
agreement is obtained between the calculated and measured diffraction intensities for this complex, 
low-symmetry surface containing only one mirror-plane symmetry element. No significant inter- 
layer spacing relaxations are found. The Debye temperature for the surface layer is found to be 
lower than in the bulk, which is indicative of larger vibrational atomic amplitudes at the surface. 
Meanwhile, the second layer shows a Debye temperature close to the bulk value. The experimental 
results for the relaxations agree well with those of our first-principles calculation. 

PACS numbers: 68.35.Bs,61.14.Hg 
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I. INTRODUCTION 



Characteristic of group V elements, bismuth crystallizes in the rhombohedral A7 structure 
as a semimetal with a small density of states at the Fermi level.— But interestingly, the 
surfaces of Bi show very different electronic properties than the bulk. Studies on the Bi(llO) 2 , 
Bi(100) 3 and Bi(lll)^ surfaces have shown that they are much more metallic than the bulk 
due to a significantly higher density of states at the Fermi level at the surface.— It has been 
found that one significant contribution is from a strong spin-orbit coupling at the surface 
due to broken inversion symmetry.— 

From a chemical point of view, the creation of a surface requires the breaking of atomic 
bonds. Covalent bonding plays only a minor role in most metals. Thus the effect of bond- 
breaking is small and surface properties are similar to those of the bulk, although localized 
electronic surface states may be present. On semiconductors, creating a surface leaves so- 
called dangling bonds which should give rise to half-filled and therefore metallic bands. 
However, it turns out that on most semiconductor surfaces the atoms re-arrange their po- 
sitions such that the dangling bonds are removed and the surface is again a semiconductor 
and not a metal.— Semimetals such as bismuth lie in between these two cases. On one hand, 
a semimetal is close to being a semiconductor since directional bonding is important and 
the valence and conduction bands are almost separated by a gap. On the other hand, there 
is a very small overlap between both bands such that the material is formally a metal. This 
delicate balance between being a metal and a semiconductor depends crucially on the atomic 
structure^ and it can be expected to be severely disturbed at the surface. 

Detailed structural information on Bi surfaces is so far limited to a recent LEED intensity- 
vs-voltage (IV) and first-principles study of the Bi(lll) surface.— One important difference 
between bulk terminated Bi(llO) and Bi(lll) is that the Bi(llO) surface exhibits dangling 
bonds, while Bi(lll) does not. In a pioneering study by Jona,— oxygen adsorption exper- 
iments suggest that Bi(llO) is noticeably more active than Bi(lll). A qualitative analysis 
of LEED patterns in Jona's study shows an unreconstructed (lxl) Bi(llO) surface struc- 
ture. From the bulk structure Jona erroneously concluded that the unit cell (and hence the 
LEED pattern) should not be exactly rectangular but that the lattice vectors should include 
an angle slightly different from 90°. This is not correct, as will become appartent below. 
The unit cell is rectangular and almost quadratic. A recent scanning tunneling microscopy 
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study by Pascual et al.fi revealed images of the Bi(llO) surface that are consistent with a 
near-square surface unit cell. 

A surface structural analysis of this surface has not been undertaken before because of 
its very low symmetry The Bi(llO) surface has only one symmetry element, a mirror plane, 
which makes the LEED-IV analysis challenging. An additional complication for Bi(llO) is 
the close stacking of atomic layers, which requires an initial calculation of the bulk diffraction 
matrices via the combined space method.— To our best knowledge, this situation has not 
been encountered in LEED-IV analyses. 

In this article, we report on a study of the surface structure of clean Bi(llO) by quantita- 
tive LEED intensity-vs-voltage analysis and ab initio calculations. Experimental diffraction 
intensities taken at a sample temperature of 110 K under normal incidence have been an- 
alyzed by comparison to dynamical LEED calculations. Great care was taken to align the 
sample considering the low symmetry diffraction pattern. The main structural parameters 
that were optimized in the LEED-IV analysis include the first 4 inter layer spacings and 
the Debye temperatures for the first 2 surface layers. Furthermore, we have performed 
first-principles calculations for the atomic structure of Bi(110). The results are in good 
agreement with the experimental relaxations. The Bi(110) surface can serve as an example 
for the application of LEED intensity analysis to low symmetry systems. 

In the following, we introduce the bulk truncated surface structure of Bi(110). Exper- 
imental and computational details are described in Sec. followed by the results and 
discussion in Sec. IIVI Conclusions are given in Sec. El 

II. BULK TRUNCATED BI(llO) SURFACE STRUCTURE 

The A7 (aArsenic) structure of bulk bismuth has a rhombohedral unit cell with a two- 
atom basis. The bulk truncated surface structure of Bi(110) is shown in Fig. ^ The side 
views show the bilayer stacking with a single bond connecting every other atom between 
neighboring bilayers. Within one bilayer, each atom in one layer bonds with two nearest- 
neighbor atoms in the other layer. The covalent bonds have been drawn by solid lines and 
the dangling bonds at the surface layer by dashed lines. The bilayer-type structure gives rise 
to alternating interlayer distances. For the truncated bulk at 110 K we have: d\ 2 = 0.208 
A, c4 = 3.064 A, d\ A = 0.208 A, d% = 3.064 A, and so on. Interlayer spacings between 
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FIG. 1: (color online) Truncated-bulk structure of Bi(llO). The solid lines and dotted lines mark 
covalent bonds and dangling bonds, respectively, (a) Top view of the first two atomic layers. Each 
layer consists of a two-dimensional rectangular lattice and the lattice constants at 110 K are given. 
The mirror planes of the structure are also shown as dashed lines, (b) and (c) Side views of the first 
eight layers perpendicular and parallel to the mirror plane, respectively. The bilayer-like structure 
with alternating short and long interlayer spacings is evident. 

the ith and jth bulk layers are indicated as g&. Noticeably, the Bi(110) surface has very 

low symmetry: the only symmetry element is a mirror plane as indicated in Fig. [ 1 The 

lengths of unit vectors at 110 K are taken as 4.731 A and 4.538 A; see Refs. ll2l)l4|l5L If 
the rhombohedral structure is treated as a pseudocubic structure as in Ref . I12I Bi(110) will 
be denoted as Bi(100). The pseudo-square character of the surface unit cell is evident: for 
a cubic Bi structure all the atoms in the first bilayer would have the same height, the unit 
cell would be rotated by about 45°, and contains only one atom. 
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III. METHODS 



A. Experimental 

The experiment was performed in a /z-metal ultrahigh vacuum chamber equipped with 
a four-grid LEED optic with a base pressure of 7xl0~ 9 Pa. Surface contamination was 
measured by Auger Electron Spectroscopy (AES) using a hemispherical electron analyzer 
and the LEED electron gun as electron source. The sample was mounted on a manipulator, 
allowing positioning to within 0.1° around all three axes of the crystal. The sample was 
cooled by liquid nitrogen. The surface was cleaned by cycles of 1 keV Ar + sputtering and 
annealing to 150°C With AES no surface contamination could be detected. The maximum 
possible oxygen contamination was determined to be 0.02 monolayers. Spot intensities were 
measured using a 16 bit Charge- Coupled Device (CCD) camera. A back-illuminated and 
Peltier cooled (— 40° C) CCD chip guaranteed a high quantum efficiency and low dark current. 
The camera was mounted on a base, which allowed rotation around all three axes. Great 
care was taken to align the camera with respect to the electron gun and the Bi crystal, as 
described below. 

To obtain intensities of the diffracted beams as a function of electron energy, the following 
procedure was employed: A series of images was recorded within the energy range from 
30 eV to 300 eV, while the energy was increased in steps of 1 eV after every recorded 
image. The integrated spot intensity of every single diffracted beam (h, k) was extracted 
from these images. The presence of only one mirror line symmetry for Bi(110) leads to 
technical challenges for the LEED experiment. These difficulties are illustrated in Fig.|2]that 
shows two measured LEED patterns taken at different incident energies. The pseudo-square 
pattern of the reciprocal lattice and the missing left/right symmetry are clearly evident. 
The up/down symmetry is given by the mirror plane in the crystal (the horizontal plane in 
Fig. I2J). It is necessary to align the sample surface perpendicular to the incoming electron 
beam, and this is usually done by comparing the IV curves of the symmetry-equivalent 
beams. Here this procedure can only be applied for the up/down angle. In order to align 
the left / right angle we optimized the diffraction spot position on the LEED screen until they 
agreed with the kinematically calculated positions. We estimate that this approach leads to 
an error of less than 1° in the angle of incidence. In the final data set, the intensities of the 
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FIG. 2: (color online) LEED patterns at two different electron beam energies for normal incidence 
on Bi(llO) at 110 K. The (1,0) and (0,1) diffraction spots are marked. 

symmetry-equivalent beams were averaged. 
B. Dynamical LEED calculations 



The dynamical LEED intensity calculations were performed using the standard package 
SATLEED (Symmetrized Automated Tensor LEED) by A. Barbieri and M.A. Van Hove^ 
within the renormalized forward scattering perturbation formalism. Atomic scattering phase 
shifts have been calculated using a muffin-tin potential model and the standard Barbieri and 
Van Hove phase shift package.— The bulk diffraction matrices for the closely spaced bilayers 
were calculated with the combined space method.^ The same muffin-tin radius of 2.87 a.u. 

n 

and phase shifts have been used as in Ref. [11|. Phase shifts have been renormalized by 
the thermal effects of root-mean-square (rms) isotropic vibrational amplitudes. Up to 15 
(L = 14) phase shifts have been used because of the strong scattering of the heavy Bi atom 
[Z = 83). The muffin-tin constant Vq is taken to be energy-independent and is optimized. 
Vi m , the imaginary part of inner potential, also referred to as damping potential or optical 
potential, is taken as 4 eV for the bulk and 4.2 eV for the first 2 overlayers. The slightly 
larger value at the surface was chosen to model the presence of dangling bonds, which 
increases the electron damping. The surface potential step of height Vq is located half a long 
bulk interlayer spacing away from the topmost layer nuclei. The bulk Debye temperature is 
fixed at 119 K,^ while the Debye temperatures for the first 2 layers are optimized. Mean- 
square atomic vibrational amplitudes < u 2 >t at temperature T for the Debye- Waller factor 
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calculation are derived from Debye temperatures G d according to the following equation:— 



where m a is the atomic mass, K Planck's constant and ks the Boltzmann constant. 

In the LEED intensity analysis, agreement between experimental and calculated LEED 
intensities is quantified by the widely used Pendry R factor, Rp, which is particularly sen- 
sitive to relative peak position and the existence of small peaks.— The uncertainties in 
the optimized structural parameters are estimated from the variation around the minimum 



where AE is the total energy range compared in the IV analysis. 19 
C. Ab initio calculations 

We have also performed ab initio calculations of the surface crystal structure of Bi(llO). 
The full-potential linearized augmented plane wave method in film-geometr y 20 ' 21 as imple- 
mented in the FLEUR-code was used and the local density approximation 22 to the density 
functional theory was employed. Spin-orbit coupling was included in the self-consistent 
calculations. 2 - The evaluation of the surface relaxation has been carried out for the sym- 
metric 14-layer film, both, with the inclusion of the spin-orbit coupling (SOC) term and 
without this term. Force calculations have been performed for the first four layers without 
spin-orbit coupling while relaxations have been carried out only for the first two interlayer 
spacings with the inclusion of SOC. In the latter evaluations we kept the interlayer spacings 
g?34 and g?45 equal to those obtained from the force calculation without SOC. The geometry 
was chosen such, that both sides of the film were terminated with an intact bilayer. A 
wavefunction cutoff of 3.8 a.u. -1 was chosen and the Brillouin zone was sampled with 32 
k-points. 
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IV. RESULTS AND DISCUSSIONS 



A. LEED structure determination 

The LEED pattern of Bi(llO) has previously been discussed by Jona.— He defined a 
pseudo-cubic bulk unit cell and concluded that the unit cell (and hence the LEED pattern) 
should not be exactly rectangular but that the lattice vectors should include an angle slightly 
different from 90°. Our study does not confirm this conclusion. Our LEED patterns as pre- 
sented in Fig. El show an exact rectangular net from careful measurements of the diffraction 
spots positions and, indeed, such an exact rectangle can also be expected from a projection 
of the bulk reciprocal lattice onto the surface.— The measured ratio of the two reciprocal unit 
cell vectors is 0.96(2) in good agreement with the expected value of 0.959. Moreover, the 
observed patterns show no indication of any reconstruction of the Bi(110) surface, despite 
the existence of active dangling bond at the surface. Apparently, Bi(110) is found to be 
very different from typical semiconductors surfaces, such as Si(100) and Ge(100) which both 
exhibit 2x1 reconstructions. 

The structural and non-structural parameters were optimized for a Bi(110) surface termi- 
nated by an intact bilayer. A termination with a split bilayer was immediately excluded due 
to lack of agreement with the experimental IV curves shown in Figs. HU and 01 20 symmetry 
in-equivalent beams with a total energy range of 3591 eV have been analyzed to determine 
the following structural and non-structural parameters: the first four interlayer spacings dij 
(j — i + 1; 1 < i < 4), the real part of inner potential Vo, and Debye temperatures 0^ 
and ©d 2 for atoms in the first and second layers, respectively. The results of the structural 
analysis are summarized in Table |U Note that the first and the third interlayer spacing cor- 
respond to the small separation (0.21 A) between the two layers making up the bilayer in the 
bulk. Their seemingly dramatic relative relaxations are very small in absolute terms. Also, 
the forth layer appears to move above the third layer by 0.01 A. However this very small 
value is clearly below our detection limit. Overall no significant relaxation for the Bi(110) 
surface is found. We have tried many possible displacement patterns allowed due to the low 
surface symmetry. However, we found no significant improvement in Rp when changing the 
relative distance between the two basis atoms in the first and second layer parallel to the 
mirror line. The Debye temperature for the first layer is found to be lower than that of the 
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TABLE I: Optimized parameter values with errors for the surface structure of Bi(llO) at 110 K 
and normal incidence. Interlayer spacings between the ith and j'th layer are indicated as dij. d\j 
is the corresponding interlayer spacing of the truncated bulk at 110 K and Ady = dij — c&. Vq 
is the real part of inner potential. @Di an d Od 2 are the Debye temperatures for the first and the 
second layer, respectively. 



Parameters Starting values Optimized values Adjj/cfe (%) 



V (eV) 


8.0 


3.5 ± 1.50 




di2 (A) 


0.208 


0.18 ±0.048 


-13 ±23 


d 2 3 (A) 


3.064 


3.06 ± 0.043 


-0.2 ±1.4 


d u (A) 


0.208 


0.01 ± 0.040 


-105 ±19 


d 45 (A) 


3.064 


3.20 ± 0.046 


±4.3 ± 1.5 


e Dl (k) 


119 






D2 (K) 


119 


H6±Sg 




Rp 




0.455 





bulk, which is consistent with an early study of Goodman and Somorjai. 2 ^ Reduced surface 
Debye temperatures are a common phenomenon reflecting the weaker bonding of surface 
atoms compared to the bulk.™ The actual numerical values of the surface Debye temper- 
ature are an important ingredient for the determination of the electron-phonon coupling 
strength from angle-resolved photoemission data . 26 i 27 i 28 Meanwhile the second layer shows 
a Debye temperature close to the bulk value. 

The LEED-IV analysis gives a relatively high Rp factor of about 0.455 compared to typ- 
ical values of 0.1 to 0.3 for clean unreconstructed metal surfaces. We believe it is due to 
the structural complexity and low symmetry of the Bi(110) surface and it is not caused by 
deviations from normal incidence during the IV measurement. We simulated non- normal 
incidence conditions extensively in the LEED-IV calculations and found that an increase in 
the incident angle gave a dramatic rise in the R factor from its minimum at zero or normal 
incidence. This shows that the sample is properly aligned. As seen on open semiconductor 
surfaces the presence of dangling bonds and the presence of voids in the open surface struc- 
ture is a real challenge for the muffin-tin approximation of the crystal potential and could 
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FIG. 3: (color online) Comparison of 20 experimental and calculated IV curves for normal incidence 
on Bi(llO) at 110 K. Solid lines show experimental data and dotted lines show calculated data. To 
be continued in Fig. 0J 



also contribute to the relatively high Rp for this surface. However, the low surface symmetry 
of Bi(110) gives rise to the large number of non-equivalent beams. The good agreement be- 
tween this large experimental data set and the calculated intensity-energy curves, as shown 
in Figs. El and 0J give us great confidence in the reliability of our results. 

The error bars of the optimized parameters were analyzed based on the variation the R 
factor around -Rp m i n , AR P = 0.043 according to Eq. |21 The dependence of Rp on a change 
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FIG. 4: Continuation of Fig. 01 

of the interlayer spacings away from their optimized values is shown in the Fig. 03 In this 
analysis, all other parameters were fixed at their optimized values. We can see that all the 
sensitivity curves take on a parabolic shape. The errors for the individual parameters are 
also listed in the Table HI 



12 



-0.05 -0.04 -0.03 -0.02 -0.01 0.00 0.01 0.02 0.03 0.04 0.05 

Deviation of interlayer spacing from best-fit value (A) 



FIG. 5: (color online) Error bar determination for the first 4 interlayer spacings based on ARp = 
0.043 and Eq. EJ 

B. Comparison to first-principles calculations 

The first-principle calculations performed for bulk Bi without the inclusion of spin-orbit 
interaction give bulk short and long interlayer spacings of 0.142 A and 3.087 A, respectively. 
Evaluations that include the SOC term lead to a very slight modification of approximately 
0.01 A of these results. Our scalar relativistic force calculations give the following values 
for the first 4 interlayer spacing relaxations at K: Adi 2 /d\ 2 = —62%, Ag^/g^ = +0.3%, 
Ad 34 /4i = -105%, and Ad 45 /d% = +4.4%. These results agree reasonably well with those 
obtained by the LEED-IV analysis at 110 K (see Table U} considering the fact that the abso- 
lute distance difference between the experimental and calculated first interlayer relaxations 
of —13% and —62% is only 0.06 A. Both the experiment and theory lead to the contraction 
of the first interlayer spacing. For the second interlayer spacing the theory gives a small 
expansion while the experiment shows a small contraction of the spacing. However the the- 
oretical result is within the experimental error bar. The absolute distance difference between 
the experimental and calculated second interlayer relaxations of 0.015 A is even smaller than 
that for the first interlayer spacing. For the third and fourth interlayer spacings the theory 
and experiment are in excellent agreement. The first-principle calculations that include the 
spin-orbit interaction term lead to Adu/d\ 2 = —43% and Ad 2 3/d 23 = +0.4% for the first 
and second interlayer spacings respectively. These values have been obtained by keeping 
the interlayer spacings Ad^/d^ and Ad^/d^ equal to those found in the scalar relativistic 
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calculations. This shows that the influence of spin-orbit interaction on the relaxation is 
small and probably will not change the values of d^ and d^ significantly. Notice, that in 
the relaxed geometry a change of 6% or 0.01 A in Ad^/d^ corresponds to an energy change 
of only 0.5 meV per surface atom which is certainly at the limit of our accuracy. 

In our force calculations, we also optimized the position of the surface atoms in a plane 
parallel to the surface. By symmetry, this movement is then confined to the mirror plane 
shown in Fig. [I] (a). We notice, that these relaxations are small and do not exceed 1.0% in 
the top four layers, consistent with the experimental findings. 

V. CONCLUSIONS 

Our results give a consistent picture of the very low-symmetry surface geometric structure 
of Bi(110) by LEED intensity analysis and first-principles calculations. Good agreement is 
reached between experimental LEED and theoretical IV curves. No structural reconstruc- 
tion occurs despite of dangling bonds present at the surface. No significant absolute value of 
relaxation is found for the first 4 interlayer spacings. The reduced top-layer Debye tempera- 
ture suggests essentially larger vibrational atomic amplitudes at the surface. Experimentally, 
the approach of sample alignment by calculating the diffraction spot positions on the LEED 
screen is very efficient and can be used for surfaces with low symmetry as well as for in- situ 
cleaved surfaces. 
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